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^ Abstract 

5^ We get deeper understanding of the role played by boundary conditions in quantum 

field theory, by studying the structure of a scalar massless quantum field theory bounded 
by two one dimensional planar crystal plates. The system can also be understood as a 
massless scalar confined to propagate in the surface of a finite cylinder. We classify the 



most general type of regular behaved boundary conditions that the quantum field can 
satisfy, in accordance with the unitarity principle of quantum field theory. Also, we 
characterize the frequency spectrum for each quantum field theory, by computing the 
holomorphic spectral function. The spectral function is the starting point to compute 
the Casimir energy as a global function over the space of allowed boundary conditions 
for the quantum field theory. 



^ 1 Introduction 
in 

^ Boundary conditions play an important role in quantum field theory. On the one hand side 
• they are generalizations of interactions which are concentrated on a region of co-dimension 1 (or 

Q higher). Typical examples are conductor boundary conditions which are a good idealization 

^ if the skin depth drops below the wavelength of the relevant electromagnetic field modes 
[T]. On the other hand, periodic boundary conditions appear naturally when considering 
topological objects, for example a quantum field on a torus, or quantum field theory at finite 

^ temperature. Other examples for the use of boundary conditions can be found in the quantum 
Hall effect [21 [3] , in the physics of graphene [H |5] , in planar plasma systems (see references 
[HI Ej) etc. Boundary phenomena are rather relevant for the structure of the vacuum and the 
low energy spectrum than for high energy phenomena. For field theories without mass scale, i. 
e. conformal theories, these effects are amplified because the long distance correlations allow 
these boundary effects to be observed over all the interior region. Precisely, the conformal field 
theories related to statistical models with second order phase transitions and string theories, 
are the ones where the boundary effects have more remarkable applications. Most of the 
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physical systems appearing in nature or implemented on a laboratory, have very large but 
finite physical dimensions. This makes of great interest the study of quantum fields confined 
in bounded domains. 

A more mathematical question is that about the most general boundary conditions com- 
patible with the general principles of self-adjointness of the Hamiltonian and unitarity. These 
questions are discussed, for example, in [8j [9j in detail. It is worth mentioning that these 
boundary conditions are in general non-local and can encode essential physical information 
to a large extend. In other words, many physical problems can be formulated in terms of 
boundary conditions. From a more mathematical point of view the boundary conditions are 
used to obtain ensure discrete spectrum of linear operators in Hilbert spaces flO\ . 

This paper is focused in the applications, of a systematic formalism for boundary conditions 
in quantum field theory [8j, to the Casimir effect. Since the seminal paper of H. B. G. Casimir 
where the effect was theoretically proposed, and the first experimental measurements in 
1958 by Sparnaay [12], many theoretical and experimental work have been done (see references 
[131 mils])- Since the theoretical and subsequent experimental discovery of the Casimir effect, 
theoretical research on Casimir effect has been focused on two separated directions. On one 
hand, the study of the dependence of the Casimir energy on the geometry of the objects 
subjected to the vacuum force. In this area are specially remarkable the references [161 HI] 
giving for the first time a general method that allows to calculate (numerically or analytically) 
the vacuum energy between compact objects of arbitrary shape like spheres and cylinders which 
is free of UV divergences. On the other hand, instead of studying the dependence of the vacuum 
energy on the geometry one may study the dependence on the boundary conditions that the 
quantum fields satisfy at the surfaces of interacting objets. Actually, related physical interest 
is the infiuence of corrugations on the Casimir force (see, for instance, in the review |14j). 
Especially periodic corrugations can be mapped onto non-trivial boundary conditions keeping 
periodicity [18]. Much work was done on periodic gratings, EH] More general methods 
developed in reference \16\ allow to implement some families of boundary conditions in the 
path integral approach. Also reference [13] shows how to introduce local boundary conditions 
on the quantum fields in the path integral approach. But the path integral approach does 
not allow to implement the most general type of boundary conditions that the quantum fields 
can satisfy (see references [2T1 122]). and especially problems when implementing Neumann 
boundary conditions are discussed in [23]). Similar problems appear in the approach [2^ 
using world-line methods. 

In the last years, the advances in experimental nano-scale physics allowed to implement a 
huge family of exotic boundary behavior, producing results of great importance in the Casimir 
effect. M. Asorey, J. M. Munoz-Castaneda et al developed in the last years a theoretical 
framework that classifies the boundary conditions that quantum fields can satisfy, as well 
as allows to study the dependence of the Casimir energy on the most general possible type 
of allowed boundary conditions (see references [H [251 ISHl |27l |2T1 [28]). The geometry of the 
bounded system in this case is a simple geometry, equivalent to the original geometry studied 
by Casimir (parallel infinite plates), but the boundary conditions are most general ones. In the 
case of infinite parallel homogeneous isotropic plates Asorey, Munoz-Castaneda et al classified 
all the allowed boundary conditions for the quantum fields attending to the type of Casimir 
force they give rise to: attractive, repulsive or null. Focused on the Casimir energy calculation, 
the central object in Asorey-Munoz-Castaneda (AMC) formalism is the spectral function. The 
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spectral function is defined as the function whose zeros are the frequencies of the quantum field 
theoretical hamiltonian of the system. Calculation of this function allowed them to compute 
the value of the Casimir energy for the most general type of boundary condition, and hence 
to study the Casimir energy as a function over the space of allowed boundary conditions. 

The purpose of this paper, is to use AMC formalism to calculate the space of boundary 
conditions in a more complicated case, and then give an explicit expression for the corre- 
sponding spectral function, in terms of the boundary condition. This calculation will provide 
the machinery necessary to compute the corresponding Casimir energy as a function over the 
space of boundary conditions that the quantum fields can satisfy. The case to be studied is 
the case of two parallel non homogenous plates (1 dimensional crystal plates), isomorphic to a 
cylindrical geometry. Via this calculation we try to get a deeper understanding of the behav- 
ior of the vacuum interaction in crystallographic systems taking into account the microscopic 
lattice structure. In order to take into account the lattice structure we will impose a periodical 
behavior over the scalar field that gives rise to the vacuum interaction. This interaction can 
be also understood as a boundary effect due to the finiteness of a 2 dimensional crystals in 
one of its directions. 

Throughout this paper we put h — c—1. 



2 The Quantum scalar field in bounded domains 

Let R X M C M.^^^ be a space-time, where the spatial surface M G M"^ is a bounded domain 
with regular boundary dM. The dynamics of a free complex scalar field propagating in M x M 
is given by the action 

5(0) = ^ f <f^\ {d^^*d,<t) + m2|0|2) - 1 / (fx^dn^, (2.1) 



where 9„ denotes the outgoing normal derivative. The boundary term is added in order to 
obtain the usual free field equations of motion in the bulk without any extra boundary term 
potential 

-(9i20 = L0, L = -A + m^ (2.2) 

where A is the Laplace operator induced in the bounded domain M by the Euclidean structure 
of R^. 

The usual canonical quantization procedure leads to the quantum Hamiltonian 

H = ^ /" (fx (t:{x)*t:{x) + 0(x)*L0(a;)) (2.3) 

where the canonical momentum density operator is given by 7r(a;) = \dt(t)*{x), and satisfies 



the canonical commutation rules 



TT X , 0[X 



-i6{x — x') 

Decomposing the fields in eigen-modes of 

L0fe(x) = Xk<t>k{x). (2.4) 

allows to write the system as an infinite collection of non-coupled harmonic oscillators. In 
order to have a unitary quantum field theory we should require the eigen-values {A^} to be 
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real (selfadjointness of L), and non negative (non-negativity of operator L). When M = M'^, 
then automatially L is self-adjoint and non-negative because A is self-adjoint and positive 
definite in M.^. But when M C M'^ is a bounded domain with regular boundary, A may not 
be essentially self- adjointj^ and in that case one needs to consider an infinite set of self-adjoint 
extensions (PIH]), which are in general not definite positive (PIE]). 

2.1 AIM formalism for self- adjoint extensions. 

To ensure self-adjointness of operator L we might use the Asorey-Ibort-Marmo (AIM) formal- 
ism for self adjoint extension^ in quantum mechanics [9]. We will study quantum mechanics 
in the domain M: 

• Quantum physical states: L^(M, C^). 

• Quantum Hamiltonian: L = — A + =^ symmetric operator in (M, C^) (smooth 

O 

functions with compact support in M), but not essentially self-adjoint in L^(M, C^). 

The obstruction for L to be self-adjoint is given by the boundary term obtained integrating 
by parts: 

(^1, L7/>2)= [ d''x4{L^2)= [ d\{L^iy^2 + t^{iJu^2). (2.5) 



S(^i, ip2) = i [(01, V2) - (v^i, V2)] dfig^ (2.6) 



being (p = ipl^^^ and = dnipl^j^^- The boundary term T,{ipi, ip2) can be used to characterize 
the domain where L is self-adjoint |9] 

• L will be self-adjoint ^{ipi, V'2) = 0. 

For simplicity reasons in this first approach, we will assume that the boundary data lie in 
L'^{dM,C). This simplification concerns also to the physical behavior of the boundary, so 
it must be carefully explained. In the most general approach cp G H~2{dM,C) and ip G 
H~2[dM, C)©ker(A"'') l^j^^, to encode all the self-adjoint extensions of A (for detailed technical 
description see chapter 1 in |8j). This restriction must be taken into account in order to allow 
non-regular behavior of the boundary data originated by physical situations that might include 
singular charge distributions on the boundary, or some kind of singular impurities distributions 
over the crystal plates. In order to simplify a first approach to the study of the vacuum 
interaction between one-dimensional crystal plates, we will avoid such kind of behaviors, and 



-'^A symmetric operator O defined on a dense domain of a Hilbert space H is said to be essentially self- 
adjoint if it admits an unique self-adjoint extension (for more details on the theory of self-adjoint operators 
and self-adjoint extensions see reference [TO]). 

^The AIM formalism for self-adjoint extensions is equivalent to the classical formalism for self-adjoint 
extensions of symmetric operators with dense domains in Hilbert spaces developed by Von Neumann (see ref. 
[lOj for modern description) in terms of defect spaces. The advantage of AIM formalism in physics is that the 
theory is formulated just in terms of boundary values (and normal derivatives) of the functions that belong 
to the domain of the self-adjoint extension. The boundary data have a clear physical interpretation, as the 
boundary data of wave functions (see [^). 
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assume for now on that all the boundary data lie in L'^{dM, C^), as was originally proposed 
in reference [9j^ Hence in this first approach we will only deal with those boundary conditions 
that do not encode physically singular behaviors of the boundaries. We will call this set of 
boundary conditions, the set of regular self-adjoint extensions (Air)- 
The first AIM theorem ensures that: 

• Self-adjoint extensions of L are in one-to-one correspondence with the unitary group 
U (L2((9M,C^)) = Mr of unitary operators over the square integrable functions over 
the boundary of the space. 

• Given a unitary operator U G A4r the associated self adjoint extension L[/ is character- 
ized by the domain of functions that satisfy the boundary condition 

ip — ii^ = U {{p + iip) (2.7) 

where now, G ^^(M, C^), and f/ G W (^^(M, C^)) 

Some examples of well known boundary conditions written in this formalism: 

• Dirichlet boundary condition: Ud = —I. 

• Neumann boundary condition: f/„ = I. 

When ±1 ^ o"(f/) unitary operator U admits non-singular Cayley transform, and boundary 
condition in this case can be rewritten in a more simpler way 

ip = Alp ; (p = A(p (2.8) 

where A is a self-adjoint operator in the domain of the corresponding self-adjoint extension. 
Cayley transformation takes unitary operators (±1 ^ cr([/)) into self-adjoint operators. 

Now we should ensure the positivity of the self-adjoint extension L^/. This condition is 
described by the Second AIM theorem, that ensures the existence of unitary operators 
f/ G for which Lu has negative energy eigen-states (also the theorem characterizes these 
operators): 

• Given U G Air with —1 G cr[U), there are near U uni-parametric families of unitary 
operators which give rise to self-adjoint extensions of L with arbitrary negative energy 
states. 

Now with these two results we should go back to the QFT and compute which of those U G Air 
give rise to self-adjoint extensions Lu that are consistent with the principles of quantum field 
theory. 

■^To include non-regular behavior of the boundary data in the formalism of quantum fields over bounded 
domains with regular boundary, and the re-formulation of the two AIM theorems see reference [S]. 
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3 Consistency conditions in QFT. 



As we saw in first two sections, a massless scalar QFT over the bounded domain M will be 
consistent if L[/ is a non negative self-adjoint extension of L. By first AIM theorem it is 
ensured that L(/ is self-adjoint for any U G M.r- It is easy to check from the expression for 
the boundary condition that: 



l-U 

l + U 



(3.1) 



By second AIM theorem there exist some unitary operators U E Ai for which (3.1) is 



negative, so not any U G Air is consistent with the principles of quantum field theory, even 
though when L(/ defines a consistent quantum mechanical system. Now observe that in right 



hand of equation (3.1): 



m 



> 



1^011 



> 



Strong Consistency Condition: we will call strongly consistent boundary conditions or 
consistent boundary conditions, those boundary conditions which give rise to a consistent 
QFT for any size 6 of the edge, i. e. those given by unitary operators U G Air satisfying 



I-U 
I + U 



^)>0, V ipEL'^idMX 



(3.2) 



That is to say that strongly consistent boundary conditions are only those characterized by 
unitary operators U G Air that satisfy 



I-U 
I + U 



> 



(3.3) 



The consistency condition can be written more explicitly in terms of the eigenvalues of the 
operator U. Since U is unitary its spectrum is a set of uni-modular complex numbers cr{U) = 
{e*^}, and taking the diagonal terms in the operator inequality given by (3.3), we get a 
condition over the eigenvalues: 



1 -e*^ 
I- > 



1 



tan 



> 



(3.4) 



Finally, with the condition (3.4) we can specify the space AirF, 



MrF = {U E M\yX = e''' e cx{U), tan ( - ) > 



(3.5) 



which is the space of consistent boundary conditions for the QFT. 
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4 2+1 QFT: 1-dimensional crystal plates. 



The analogue of the 1+1 dimensional example explained above is the case of a 2+1 dimensional 
quantum field theory, in a bounded space limited by two infinite homogeneous isotropic parallel 
wires. This case and its generalizations to general D+1 scalar quantum field theories has been 
studied in detail in [8] . The simplest extension that can be done for all these plate geometries 
is to broke the isotropy of the plates, by introducing a D-1 dimensional lattice over each plate. 
To learn how the boundary conditions behave in this case the simplest case we can study is 
that one given by two parallel and identical crystal wires, as is shown in figure [TJ 
We are going to study the system of one scalar field governed by the action 



S{<P) 



'(9„0 + 



KxSM 



d X(f)*dn(, 



(4.1) 



with physical space given in figure [Tj From geometrical point of view the physical space where 
the quantum field is confined is M = [0, L] x M, but there are two identical lattices at x = 
and X = L respectively. This fact implies that, the boundary data of the quantum field should 
preserve this lattice structure, i. e. the boundary data should have the same periodicity of the 
lattice. Since the two lattices are identical, and no displacement between them appears, in this 
case boundary data should only have an unique periodicity. This geometry is isomorphic to a 
cylindrical space-time where the longitudinal direction corresponds to the direction orthogonal 
to the plates, and the circled direction is the parallel direction to the plates. In addition, the 
simple periodic structure imposed implies that the field over the cylinder has winding number 
1. 



x = 



M 



X = L 



Figure 1: The physical space hmited by two one dimensional lattices. 



In order to determine the set AirF of regular self-adjoint extensions, we need to construct 
the most general wave function for the quantum mechanical system associated to the QFT 
between crystal plates. Since there is no potential between plates, and in the plates, the only 
requirement to be imposed, is the periodicity of the wave function the direction parallel to the 
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crystal (the y-direction) : 



?/'„,fe(x, y) = Ae'('=^+^"^) + BS^^'-'i-y^ + Ce^(-^^+'?"^) + De-^('^^+«"^) (4.2) 

being g„ = 27m/a with n G N. Since the momenta in the y-direction are discrete (due to 
the requirement of periodicity imposed in order to take into account the crystal structure 
and Bloch's theorem), and the x direction is compact and bounded (x G [0,L]), the allowed 
momenta in the x-direction are also a discrete set. The allowed values of momentum k should 
arise from the boundary condition. 



5 The set of regular boundary conditions: Aij.p. 



In order to exactly determine the set Air of regular boundary conditions (reduced self-adjoint 
extensions), we need to explicitly compute the equations that arise from the general boundary 
condition given by equation 2.7 Using same notation used in [Hj the boundary data can be 
written as 

V'fcn(0,i/) =F idxi>kniO,y) 
and the boundary condition takes the form 



(5.1) 



(5.2) 



From equation 4.2 we can obtain an explicit expression for \E' 



■ipk,n{L, y) ± idxiJk,n{L, y) 



+ {B{l±tk) + D{lTik))e-'"'-y 
= {Ae'''^il±ik) + Ce-"'^ilTik))e 
+ {Be'''^{l±ik) + De-'''^{lT^k))e- 



1'gny 



(5.3) 



(5.4) 



In view of these last expressions, and taking into account that the functions e*''"^ and e"*"^"^, 
are linear independent, we can rewrite all in terms of four dimensional column vectors and 
4x4 matrice^ The 2 dimensional column vector ^>± becomes now 



A{l±ik) + C{lTik) 
B{l±ik) + D{lTik) 
Ae'''^{l ± ik) + Ce-'''^{1 t ik) 
\ Be^^^ll ± ik) + De-'^^{1 t ik) J 



(5.5) 



and the boundary condition is written in terms of a 4x4 unitary matrix: 

^.{k-L)- = u(^^{k-L))- f/GU(4) 



(5.6) 



^{^±^ are the e*''"'^ components, and (^^±) the e components. Also we could have imposed 

the convention of choosing (^^±) to be the e*'"^ components, and (^^±) the e"'*"^ components, but 

the spectral function does not change under this change, due to the properties of the determinant (det(j4B) = 
det(^) det(B)). 
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Observation. Not all matrices U G U(4) are allowed boundary conditions. This means, that 
not for all U G U(4) the theory of quantum fields obtained is unitary, and well defined. 

In order to classify whicli U G U(4) give rise to unitary QFT's we must make use of the 
consistency lemma (see [8]). From the consistency lemma is rapidly deduced that the space 
of regular boundary conditions is given by the subset: 



Mr = |f/ G U(4)| Ve*'' G spec([/),tan > Oj (5.7) 

A^r- is a space with boundary, given by: 

dMr = |f/ G U(4)| 3e'^ G spec(f/),tan = o| (5.8) 

Hence, Air is a 15 dimensional manifold, with 14 dimensional boundary dAir- Since the 
sequence 

W U(l) U(4) ^ U(l) ^ 1 (5.9) 

is exact via the determinant projection of U(4) into U(l), we can write U(4) = U(l) x SU(4) 
to get a simphfied coordinate system in U(4), and hence also for M-r- 

U(4) = {e'^U\U G SU(4), e*^ G U(l)} (5.10) 

An explicit parametrization for SU(4) will be given later. 



6 Explicit computation of the spectral function. 

In this section we will compute explicitly the spectral function whose zeroes give the allowed 
momenta for the transverse direction, in terms of the boundary condition. This computation 
has never been done before, and is very useful from the computational point of view. 
Defining the four-vector 

/ A \ 

B 



C 



and having in mind equation (5.5), rapidly allows us to write as 

*± = M±(A;)<i>. 



The matrices M±{k) are defined as: 

/ 



l^ik 




iak I 



\ 



l±ik) 






l^ik 


3^'^^(1 ±ifc) 



l±ik 








l±ik 




e-^a>'(^lzpik) j 



Now the boundary condition is written in terms of M±{}z) and $ as 

[M_(A;)-[/-M+(A;)]-<l> = 



(6.i; 



(6.2) 



(6.3) 



(6.4) 
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Observation. The matrices M±{k) have determinant 

D± (k) = -4 [2k cos{kL) ^ - 1) sm(A;L)] ^ (6.5) 

and therefore they are singular for a discrete but infinite number of values of k, i. e. the 
equation D±{k) = has an infinite number of solutions. Nevertheless, D±{k) = does not 
determine the allowed values of k because roots of D^{k) do not coincide with roots of D_{k). 



In order to have non-trivial solutions satisfying (6.4), and in view of the observation made 
above, we must require the condition: 

det [M_{k) -U ■ M+{k)] = (6.6) 

This condition, as is known from |8j gives rise to the spectral function, that is now defined as: 

hu{k; L) = det [M_{k) - U ■ M+{k)] (6.7) 

Lemma 1. The infinite discrete set of allowed values for the transverse momentum k are 
given by the set of real positive zeros of the spectral function hu{k; L) . Note that AIM theorem 
ensures the self-adjointness, meanwhile the consistency lemma ensures the unitarity of the 
QFT. This means that after imposing both conditions, the boundary condition can not give 
rise to bound states, i. e., these two conditions ensure that the spectral function does not have 
imaginary zeros. 

Explicit computation of the spectral function can be very useful form a pure computational 
point of view. Even when the spectral function has a very large expression, is fundamentally 
composed by trigonometric functions with argument kL and polynomials in k whose coeffi- 
cients depend on the boundary condition. This means that calculating its zeros numerically is 
not a difficult task, with most recent methods. We proceed now to give the explicit algebraic 
structure of the spectral function. Later on some particular cases for matriz U will be studied. 

For studying the algebraic structure of hu{k;L) we introduce an auxiliary real positive 
parameter e, and define 

hu{k; L- e) == det [M_(fc) - eU ■ M+{k)] . (6.8) 

Trivially we can re-store the spectral function by making e = 1. The key point is that 
hu{k;L]e) is a four order polynomial in e, and the order in e coincides with the number of 
[/-matrix elements appearing. Hence we can rapidly identify where the algebraic invariants of 
matrix U appear: 

• Terms of order 4 in e are the only ones that include det(f/). 

• Terms of order 1 in e are the only ones where tr(f/) can appear. 

• Terms of order in e do not depend on matrix U and hence are independent of the 
boundary condition. 

We introduce at this point the notation for the coefficients in e of hu{k; L; e): 

hu{k; L; e) = co{k, L) + ci(A;, L, U)e + C2(A;, L, U)e^ + c^{k, L, U)€^ + c^^k, L, U)e^ (6.9) 

Coefficients itemized before as the ones that include algebraic invariants of the U matrix, will 
be named from now on, algebraic terms of hu{k; L; e) . 
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6.1 The algebraic terms of hu{k;L;e). 

In this sub-section we show exphcitly how the the algebraic invariants of matrix U enter in 
the spectral function via the hu function. 



The coefHcient Co(A;, L) intuitively can only depend on the matrix M_ since it is the only 
part of hu{k; L; e) that does not multiplied by U. The explicit calculation ensures that intuition 
points in the right direction, and the order coefficient is given by: 

co{k,L) = det (M_(/c,L)) = -4 [2A;cos(A;L) + {k^ - l)sin(A;L)]^ (6.10) 

The coefficient ci{k, L, U) is the one where tr(C/) can appear. In general also non-diagonal 
terms can appear. After calculation, the final expression for the coefficient ci(A;, L, U) is: 

ci{k,L,U) = {2kcos{kL) + {k^-l)sm(kL))x 

+ l)tr([/) sin(A;L) + 8k {U^ + U31 + U24 + (6.11) 

As can be rapidly seen, ci{k,L,U) depends on tr{U) but also has a non-algebraic part de- 
pending on non-diagonal matrix elements. 



The coefficient C4{k, L, U) is the one where the determinant appears. Since it comes from 
the fourth order in the e expansin, in this case no non algebraic terms can appear, since all 
combinations of four matrix elements of U appear in det([/). After calculation and simplifi- 
cation, the corresponding coefficient reads: 

C4(A;, L, U) = det(C/)dei(M+) = -4 det(C/) [2A; cos(A;L) - (A;^ - 1) sin(A;L)] ^ (6.12) 
6.2 The non-algebraic terms of hu{k;L;e). 

Calculation of the non-algebraic terms is a more tedious than for the algebraic ones. Never- 
theless they can also be computed using standard simplification techniques. 



The second order coefficient C2{k,L,U) contains only combinations of two matrix ele- 
ments of matrix U. We could in general expect all possible combinations to appear. After 
calculation and simplification the coefficient is given by: 



C2{k, L, U) = we (f/i4f/23 - U13U24 - U21USI + U21UM + Us2Uu + UM2 - UM2 + Ui2Ua^) 

+ 4 + l)\m\kL) {UMi - U^iU22 + U23U32 - ^^22^33 

+ UuU^i + Us^U^s + UnUu - U33U44) 

+ (l6k^ - 4 (A;2 + l)\m\kL)'^ {U^^U^^ - + U22UU - U24U42) 

- 8k {e + 1) sin(A;L) (C/12C/23 + C/34C/23 - C/iiC/24 - C/22C/31 

+ U21U32 - U24U33 + U12U41 + U34U41 - C/iiC/42 - U33U42 + U32U43 + U14U21 

+ C/14C/43 - C/31C/44 - C/13C/22 + C/13C/44) (6.13) 
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The coefRcient C3{k, L, U) is obtained after a similar calculation of the one made for C2'- 

csik, L, U) = -4(A';2 + 1) sin(A:L) ((A;^ - 1) sin(A;L) + 2A; cos(A:L)) [UUU22U33 

+ U12U23U31 + U13U32U21 — U11U32U23 — U13U22U31 — U12U21U33 

+ U11U22U44 + U12U24U41 + U21U42U14 — U11U24U42 — U14U22U41 

- Ui2U2iU4i + UnUssUu + U13U34U41 + U31U43U14 - U11U34U43 

- U14U33U41 - U13U31U44 + U22U33U44 + U23U34U42 + U32U43U24 

- U22U34U43 — U24U33U42 — U23U32U44] 

- 8k {2k cos{kL) - (A;^ - 1) sin(A;L)) [UUU33U24 + U11U33U42 

+ U13U32U41 + U13U34U21 + U31U12U43 + U31U14U23 - U11U23U34 

- U12U33U41 — U11U32U43 — U21U33U14 — U13U31U42 — U13U31U24 
+ U22U44U13 + U22U44U31 + U24U41U32 + U42U21U34 + U24U43U12 
+ U42U23U14 — U22U14U43 — U22U34U41 — U12U44U23 — U32U44U21 

- U24U42U13 - U24U42U31] (6.14) 

Although C3{k, L, U) can be though to be constructed in terms of non algebraic invariants of 
matrix [/, it can be rewritten in terms of algebraic invariants of a matrix associated to the 
matrix U : the matrix of minors. Also this way of writing C3 provides a way to simplify a 
lot the expression for C3. First of all recall the definition of the matrix of minors. In the most 
general case, if f/ is a square nxn matrix the n'^ (n — 1) x (n — 1) submatrices are constructed 
deleting one row and one column, and we use the notation: 

jjihj) ^ {delete row n + 1 — i, and column n + 1 — j} (6.15) 

With this notation, we define the nxn matrix of minors A^^^ associated to U, as the matrix 
whose elements a^^^ are given by: 

ajf = det (C/(^'^)) (6.16) 

Using the matrix of minors associated to U, the coefficient C3 can be written in a very compact 
and symmetric way: 

C3(A;, L, U) = {2k cos{kL) - {k^ - 1) sin(A;L)) [4{k^ + 1) sin(/tL)tr {A^^^) 

+ 8;t(a!? + 4r) + 4? + al?)] (6.17) 

After this simplification, we have explicitly computed the expression for the spectral func- 
tion whose zeros are the allowed values for the transverse momentum: 

hu{k; L) = co(A;, L) + C4 {k, L, det{U)) + C3 {k, L, A^^^) + C2{k, L, U) (6.18) 

7 Particular cases. 

7.1 Diagonal, and anti-diagonal elements of M.r' 

Diagonal elements of Mr are given by four parameters, that are their eigenvalues: 

C/d = diag(Ai,...,A4) (7.1) 
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Each eigenvalue is a modulus 1 complex number Xj = e^^ , satisfying the consistency condition 
tan {9j/2) > 0. In this case the corresponding minors matrix A^^''-^ is also diagonal, and given 
by: 

>1(^<^) = diag (Ai A2A3, A1A2A4, A1A3A4, A2A3A4) (7.2) 
If we define the two variable polynomial 

p{x, y; k, L) = 2i {k'^ + l) sin(A:L)(a; + y) - {D+f'^xy + {D^f'^ (7.3) 

the corresponding spectral function can be written in a compact form as: 

hui{k, L) = p(Ai, A3; k, L)p{X2, A4; k, L) (7.4) 

There are some particular cases of diagonal boundary conditions that should be explicitly 
computed. The two most important cases are Dirichlet and Neumann boundary conditions, 
given by: 

Ud = -I; Un = I (7.5) 

The spectral function for both cases looks very simple, and its roots can be calculated analyt- 
ically: 

hoik^L) = -%Asm^{kL), hoik^L) = -Uk^ sm^{kL) (7.6) 

Also we can mix Neumann and Dirichlet boundary conditions, by imposing one boundary 
condition for the upwards modes (e*^^ components), and the other one for the downards 
modes (e~'*^ components), which is given by a diagonal matrix 

u.o ( ; ) (7.7) 

In this case the spectral function resulting is also very simple, and given by 

hND{k, L) = -64A;2 cos^{ak) (7.8) 



Anti-diagonal elements of Mr can be handled with analytical techniques. Let us write 
an arbitrary anti-diagonal element of A4r as 

/ Ai \ 
A2 

A3 

\ A4 y 



u, 



ad 



The corresponding matrix of minors is easily computed, and has the form: 

/ 










y — A1A2A3 






- Ai A2A4 






-A1A3A4 





"A2A3A4 \ 






(7.9) 



(7.10) 



Simplifying the general expression for the spectral function in the case of anti-diagonal bound- 
ary conditions, we obtain the spectral function 



had{k, L)^D_ + L>+AiA2A3A4 + A{e + if sin\kL){X2X3 + A1A4) + 16A;'(AiA2 + A3A4) (7.11) 
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Well known boundary conditions belong to this group. The most important ones are periodic, 
and anti-periodic boundary conditions. For the periodic boundary conditions the correspond- 
ing matrix is 



( 








1 \ 








1 








1 








^ 1 











(7.12) 



and the spectral function obtained has the form: 

hp{k, L) = 64A;2 sm^{kL) (7.13) 

For the case of anti-periodic boundary conditions, the corresponding matrix is Uap — —Up, 
which leads to the same spectral function as for periodic boundary conditions: hap{k, L) — 
hp{k, L). We can also mix periodic and anti-periodic boundary conditions making the system 
to distinguish between upwards propagating modes, and downwards propagating modes: 



/ 











1 \ 










1 










-1 








V 


-1 











(7.14) 



The spectral function generated in this case is given by 

hpapik, L) = -16(-1 + k^ sin^ikL). (7.15) 

For this case is immediately observed from the spectral function that there appears a doubly 
degenerated state for k = 1 which does not depend on the distance L between crystal wires, 
and hence does not contribute to the vacuum force between plates. 

7.2 Box-anti-diagonal elements of A4r. 

Other well known important boundary conditions, studied for the case of homogeneous isotropic 
plates, arise in the case of one dimensional perfect lattices, as box-anti-diagonal boundary con- 
ditions, i.e. boundary conditions given by matrices of the form: 

where A,B E U(2). This subset of boundary conditions is contained in (U(2) x U(2)) P) M.rF- 
The two most well-known boundary conditions belonging to this subset arc the pseudo-periodic 
and quasi-periodic boundary conditions. Both examples correspond to one-parameter set of 
boundary conditions that interpolate between periodic and anti-periodic boundary conditions. 

The pseudo-periodic boundeiry conditions physically correspond to a situation in which 
both plates are identified (cylinder geometry), with a magnetic flux passing along the cylinder. 
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This situation can be developed with a boundary condition given by the box diagonal matrix 



/ 











e 

ia 



■la 



e"*" 
e^" 

0/ 



(7.17) 



It must be pointed, that if the upwards modes are suffering a flux a then, the downwards modes 
coupled to the same flux are suffering a flux —a, as is perfectly shown in the corresponding 
unitary matrix. The resulting spectral function is given by 



hpp{k, L) = 8(-l + Qk^ - r + (1 + ky cos(2a)) sin^(A;L) 



(7.18) 



If now we couple downwards and upwards modes to different fluxes, the corresponding bound- 
ary condition is given by 



mpp 



( 





i/3 










e"*° 
e*" 

0/ 



(7.19) 



which leads to the spectral function 

h^ppik, L) = 8(-l + 6k^ -k^ + {l + k^ cos(a + /?)) sm'^{kL) 



(7.20) 



As expected, when /3 = a the periodic boundary conditions are re-storcd, because the flux 
suffered by upwards modes, is compensated by the effect of the flux suffered by the downwards 
modes. 



Quasi-periodic boundairy conditions are mathematically very similar to pseudo-periodic 
boundary conditions, but the corresponding physical situation is very different. Quasi-periodic 
boundary conditions correspond to an identification of both plates via a delta function, i.e an 
step discontinuity in the first derivatives at the junction point. These boundary conditions 
are given by the matrix 



U, 



qp 






cos(a;) 
y — sin(Q; 

and the corresponding spectral function reads: 






— sin(Q;) 

— cos(a) 



cosya ) 
sin(Q;) 







sml^aj 
— cos(q;) 





(7.21) 



hqp{k, L) = 8(-l + 6r - r + (1 + ky cos(2q;) sin^(A;L) 



(7.22) 



8 Conclusions. 

In the first part of the paper we classified the set boundary conditions that give rise to a unitary 
massless scalar quantum field theory with regular behavior of the field in the boundary (given 
by the two one dimensional crystal plates, or equivalently by two circles). In the second part 
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we determined uniquely the spectrum of L„ by calculating the spectral function for each self 
adjoint extensin Lijj G AirF- The importance of the spectral function in relation to Casimir 
effect is that having the spectral function explicitly computed it can be computed the Casimir 
energy as a global function over the 16 dimensional space AirF- Taking into account that the 
vacuum energy is given basically by 



a2 



and the fact that the spectral function hu{k) is holomorphic in k, we can replace the summation 
in k by an integration in the complex /c-plane over a contour C that encloses all the zeroes of 
hu{k): 



kez{hu) 

Hence using the heat kernel regularization 



EAU^,.)^J: ^ J,. + l!j!e-V^^ (8.3) 

in the same way as done in [8] the spectral function allows to obtain a computable (numerically 
in general and analytically in particular simple cases) expression for the vacuum energy and 
the corresponding force. 

Finally in the third part of the paper, we have computed explicitly the spectral function 
for the most well-know types of boundary conditions, as well as some generalizations to one- 
parameter and two-parameter families of boundary conditions. 
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